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Abstract 

Asymptotic states in long Josephson junctions are investigated in 
an external magnetic field. We show that a choice one of the solution 
of the stationary Ferrell-Prange equation can carry be out with use of 
an asymptotic solution of the sine-Gordon equation and that an evo- 
lution to that stable solution occurs by passing through metastable 
states, which is determined with a form of quickly damped initial per- 
turbation. The boundary sine-Gordon and Ferrell-Prange problems 
were carried out with a numerical simulation. An approximated ex- 
pression for the vortex and antivortex states is obtained in the case of 
large values of an external magnetic field. 

The problem of magnetic field penetration in long Josephson junctions is 
in a sense a classical one Nevertheless, interest in it has not abated 

[P]-|13|. For example, in long Josephson junctions a dynamical chaos was 



discovered p|-||TO|. A long Josephson junction is a very good system to 
study Josephson vortex (fluxons, solitons) motions and vortex interactions 

m-m- 

It is well known that a weak magnetic field penetrates exponentially in a 
long Josephson junction. In a strong magnetic field a state of the long Joseph- 
son junction arises that is characterized by the appearance of a Josephson 
vortex. An equation that describes this phenomenon is nonlinear and in a 
general case it is the driven damped sine-Gordon equation. In a station- 
ary case we have the nonlinear Ferrell-Prange equation and one would think 
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the stationary states of a long Josephson junction represent solutions of the 
above equation, although it turned out that the Ferrell-Prange equation has 
more than one solution and the number of them increases by increasing the 
total length L of the long Josephson junction. Consequently, the problem 
of the selection of solutions arises. Considering this problem from the ther- 
modynamic point of view one can affirm that a solution is realized which 
corresponds to a minimum of the thermodynamic potential. However, the 
situation turns out to be more complicated. In fact, if we suppose that the 
minimum of the thermodynamic potential is satisfied with some solutions 
of the stationary equation, then this criterion for selection of the solutions 
turns out to be insufficient. Apparently this is precisely what takes place in 
our case. In the present work it is shown that the problem of selection of so- 
lutions for the stationary Ferrell-Prange equation is unequivocally solved in 
the following way: An asymptotic solution of the nonstationary sine-Gordon 
equation is found that coincides with the solution of the stationary Ferrell- 
Prange equation by t ^ oo, and what is more either of asimptotic solution, 
which is realized, depends on a form of rapid damped initial perturbation 
or on a way of switching on an external magnetic field. This exceptionally 
important role of initial perturbation that appears as a trigger mechanism is 
obviously explained by system nonlinearity. 

We consider a long Josephson junction in an external magnetic field Hext 
perpendicular to the junction. The Josephson phase variable (p{x,t) in a 
long Josephson junction is satisfied with the one-dimensional sine-Gordon 
equation 



where x is the distance along the junction normalized to the Josephson pen- 
etration length Xj, 



t is the time normalized to the inverse of the Josephson plasma frequency 



7 is the dissipative coefficient per unit area, $0 is the flux quantum, jc is 
the ctitical current density of superconductor, d = 2XL + b, Xl is the London 



(fitt{x, t) + 2j(ft{x, t) - ipxx{x, t) 



sin ip{x, t), 



(1) 
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penetration length, 6 is the thickness of the dielectric barrier, C is the junction 
capacitance per unit area. 

Conditions corresponding to the presence of the external magnetic field 
can be modeled by the following boundary conditions: 



Here and below the magnetic field is normalized to $o/27rAj(i, and the total 
length of the junction L is normalized to Aj. 

First of all we consider an asimptotic solution (by t ^ oo) oi Eq. (1) 
with the boundary conditions (2) and suppose that limt^^ Hgxt{0,t) — Hq. 
At first sight this problem is trivial, since all asymptotic solutions will be 
stationary on account of dissipation, they will be satisfied with the stationary 
Ferrell-Prange equation with the corresponding boundary conditions: 



However, numerical simulations of the boundary problem (1),(2) and (3) 
show that not all solutions of the problem (3) are asymptotic solutions of 
the problem (1), (2) or in other words we can say that not all solutions of 
the problem (3) are physically observed. For example, for the small external 
fields Hq < 1 there exists a solution of the problem (3) that corresponds to 
an observed dumping of the magnetic field into a junction. Simultaneously 
with that solution another solution of the problem (3) exists that gives very 
strong growth of a field in the same direction into a long Josephson junction. 
To cut off "nonphysical" solutions supplementry considerations arc attracted 
in the form of a condition of the minimum of a thermodynamic potential, as 
was mentioned above. 

The results of problem (1), (2) simulations showed that "nonphysical" 
solutions are never realized in the asymptotic approximation. The reason for 
the lack of the asymptotic "nonphysical" solutions is their instability relative 
to small perturbations. 

We show that the stationary solution ipo{x) of the problem (3) is stable 
then and only then when all values of the spectrum E are positive, i.e., 
-E'min > 0, for the boundary problem 



(2) 



(3) 



Uxxi^) + u{x) COS (fio{x) — Eu{x), 



(4) 



Ux{x)\x=0 = Ux{x)\x=L 



0. 
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To prove this assertion we shall linearize the sine-Gordon equation in 
the vicinity of stationary solution (Pq{x); i.e., we shall suppose (p{x,t) = 
^o{x) + t), where 9{x, t) is an infinitesimal perturbation. We obtain the 
equation for the function 9{x,t) from the sine-Gordon equation taking into 
account the Ferrell-Prange equation: 

9tt{x, t) + 2-f9t{x, t) - 9^^{x, t) = -9{x, t) cos (po{x), 

(5) 

9x{x)\^=o = 9x{x)\x=L = 0. 

We can obtain a general solution of the boundary Eq. (5) by means of the 
expansion of the function 6{x,t) in a series in terms of a complete system of 
eigenfunctions of the Schrodinger operator with the potential cos<^o(2^): 

^(x,i) = 5^e^"V(^), (6) 

n 

where A„ and Un{x) are eigenvalues and eigenfunctions of the Schrodinger 
operator of the problem (4), respectively. Substituting expansion (6) into 
Eq. (5) and taking into account Eq. (4), we get for A^: 

A„ = -7±^72-£;„. (7) 

It is clear that if even one of the numbers A„ has a positive real part, then the 
infinitesimal perturbations 9{x, t) will increase exponentially in the course of 
time and, further, if all A„ have a negative real part, then the perturbations 
damp exponentially. Since the Schrodinger operator is a Hcrmitian operator 
in a space of functions which have a derivative on the ends of intervals equal 
to zero, its eigenvalues En are real. It results from expression (7) that A„ 
can have a positive real part only with En < 0. Thus with En > all A^ < 
and perturbations 9{x, t) damp exponentially. It is clear that the condition 
Enmin > is the condition of the stable solution of the problem (3) with any 
values 7 > 0. 

One can call unstable solutions (po{x) of the problem (3) metastable, 
since they may be long-lived. Analysis of the spectrum problem gives in- 
formation about the survival time of metastable states. Let 3?A„ > 0; then 
T = {^Xn)~^ makes sense of the characteristic decay time of the metastable 
state. The corresponding eigenfunctions Un{x) are perturbations that destroy 
the metastable states. 
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A detailed analysis of the problem (1), (2), (3), and (4) is failed to be 
carried out by analytical methods. The numerical simulation results of this 
problems showed the following. 

(1) The number N of stable and metastable states, i.e., the number of 
solutions of the problem (3), increases with increasing L. For example, with 
i/o > 2 for L = 27r, = 4, for L = 20, N = 8. 

(2) From all solutions of the problem (3) approximately half of the states 
is stable. For example, with Hq — 3 ior L — 27r, 2 states out of 4 are stable; 
for L = 20, 3 states out of 8 arc stable; for L = 40, 11 states of 18 are stable. 

(3) The numerical investigations of the sine-Gordon Eq. (1) showed that 
either stable state [stable solution of the problem (3)] is realized depending 
on the method of superposition of an external field in a long Josephson 
junction. In other words, an established magnetic field and currents in the 
junction "remember" how an external magnetic field was started. 

For the case Hq = 2, L = 27r switching on the magnetic field in the form 
-f^ext(0,t) = HcxtiL,t) = Ho[l — aexp{—t/5)cost] leads to the system coming 
to one or another state (Fig. 1) depending on the controlling parameter a. 
So with < a < 0.36 state 4 is reahzed in Fig. 1; with 0.37 < a < 0.74 the 
one fluxon state 6 is realized; with a > 0.75 the two fluxon state 8 is realized. 

It turns out that transition to one of the stable asymptotic states occurs 
through the metastable states. For example, with a = 0.74 state 6 arises 
from the metastable state 7; with a = 0.75 state 8 arises from the same 
metastable state 7 (Fig. 1). This property of the system is not characteristic 
for the case Hq — 2, corresponding to a separatrix. For other values of Hq 
the same property takes place. In Fig. 2 the evolution of the magnetic field 
is shown for the case Hq = 2 with different values of the parameter a into 
the junction. 

In our opinion the conclusion is the most important. Thus, an asymp- 
totic state of a magnetic field, realized in a long Josephson junction depends 
essentially on a method of switching on an external field, which is a direct 
result of the sine-Gordon equation. The asymptote of problem (1), (2) taking 
into account the method of switching on an external field does select one of 
the stable solutions of the stationary Ferrell-Prange equation. It is possible 
that this property of the solution is general enough for nonlinear systems. 
With the condition of minimum thermodynamic potential all stable solutions 
of problem (3) are satisfied apparently. 

We should like to note another interesting fact following from the nu- 
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merical simulation of problems (1), (2) and (3). Among all the stationary 
solutions there exists one which corresponds to the maximum of a magnetic 
field and the another that corresponds to the minimum of a magnetic field. 
In Fig. 1 these are states 8 and 6, respectively. The states are stable. State 
8 in this figure is called a vortex; state 6 in the same figure is called an 
antivortex. For these states wc can obtain an approximate formula solving 
problem (3) by means of the itcraction method 

Voli^) ^s^^Vo'^\x), (8) 

^Ox \x=0 — ^Ox \x=L — -no- 
Suppose ipQ\x) — Hqx + a, where a is some kind of constant. Then 

Voxxix) = sm{Hox + a), 



Vox \x=0 — <fOx \x=L — -Ho- 



(9) 



Integrating (8), we find 



= v9S(x) = Ho±^ fcos(i/oa; - ^) - cos(:^)l . (10) 



H, 



2 ' '2 



In formula (10) the "+" sign corresponds to a vortex state; the "— " sign 
corresponds to an antivortex state. Expression (10) is correct in the case of 
large fields. Comparision of the values, obtained by means of formula (10), 
with results of numerical simulation shows that they practically do not differ 
for = 10. 

In conclusion we note that a swing of deflections of the field iJ*^^-* [x] from 
Hq makes 2/ Hq and this swing decrease by increasing of the external field. 
In other words one can say that the number of the stationary states per unit 
of field increases. In turn this circumstance may lead to the small fluctuations 
of the external field Hq to cause of a "jump" from the one stable state to the 
other. We hope to discuss this problem in detail in our next work. 
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Fig. 1. Solutions of boundary problem (3) for the case Hq = 2, L = 27r. Solid 
line (4,6,8), stable states; dashed line (1,2,3,5,7), metastable states. 

Fig. 2. Evolution of magnetic field for the case Hq = 2, L = 2ti with different 
values of parameter a. (a) a = 0.2, (b) a = 0.5, and (c) a = 1.0. 



